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An explicit fluid flow simulation of electromagnetic wave propagation in the gravita- 
tional field of a Schwarzschild black hole is given. The fluid has a constant refractive 
index and a spherically symmetric inward directed flow. The resulting form of the 
\ metric leads to a new coordinate system in which the Schwarzschild vacuum is written 

in Gordon's form. It is shown that a closely related coordinate system interpolates 
■^j- ' between the Kerr-Schild and Painleve-Gullstrand coordinates. 
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1 Introduction 

Analogue models of gravitational fields have been rather extensively discussed in recent years (see e.g. £Q). 
One motivation is that it might become possible in the near future to simulate kinematic aspects of strong 
gravitational fields. A number of different ways of achieving this goal have been proposed, mainly by using 
acoustic or electromagnetic waves. An effect which could in principle be measured by simulating a black hole 
would be the Hawking radiation associated with the black hole horizon [2] . An encouraging development is the 
recent experiments on media with a very large effective index of refraction 3 . Since the light speeds in these 
experiments are of the order of meters or centimeters per second, it would be possible to simulate an horizon 
which could lead to observable Hawking radiation. The horizon will be the boundary of the region where the 
fluid velocity is larger than the velocity of light in the medium. A complication is that the low values of the 
light velocity are due to a large dispersion. Therefore, it will not be possible in general to give a geometric 
formulation of the dynamics of light propagation in such a situation. However, if the dispersion is sufficiently 
"well-behaved", it seems that one can still define an effective metric for light propagation 0]. In this note we 
will consider only the more simple dispersion-free case. 

We will present an explicit theoretical model for simulating the Schwarzschild vacuum field, using a (hy- 
pothetical) dispersion-free transparent dielectric fluid. The refractive index is assumed to be independent of 
the density. The fluid flow is spherically symmetric and inward directed. A byproduct of this work is a new 
representation of the Schwarzschild gravitational field in terms of the Gordon form corresponding to this flow. 
The metric form given by Gordon is the natural representation of the (simulated) gravitational field in 
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this context. It is this particular way of writing the metric which results in a new coordinate system for the 
Schwarzschild vacuum. Like some other coordinate systems (e.g. Kruskal @, Painleve-Gullstrand |S] and 
Kerr-Schild [Hj) it is non-singular at the horizon and can be continued all the way down to the central singu- 
larity. Another nice property of this representation of the gravitational field is that it can be viewed as a finite 
perturbation of the Minkowski space-time. It shares this property with the Kerr-Schild form of the metric. 
However, it has the additional advantage of having the direct physical interpretation in terms of the flowing 
dielectric fluid. 

Consider now a transparent fluid moving with 4- velocity in our laboratory. In practical applications, 
the velocity of the fluid will always be non-relativistic. It is nevertheless useful to use a relativistic notation in 
this context. The lab geometry is given by the usual Minkowski metric, here shown in Lorentzian coordinates 
x» = (ct,x,y,z) 

r) = -c 2 dt 2 + dx 2 + dy 2 + dz 2 . (1) 

If the fluid has a refractive index, n, which is constant or slowly varying (in space and time), then the propagation 
of electromagnetic waves is governed by a wave operator 

d 2 

9 " dx^dx" ' ^ 

where g^ u is Gordon's metric [5] 

/^f -m«V , (3) 
with m = n 2 — 1. The covariant form of Gordon's metric is 



SV = V»v + k , (4) 



where k = m/n = 1 — n and = rj^u 



u 



2 Simulating the Schwarzschild vacuum geometry by a moving fluid 

We will now show how the Schwarzschild vacuum can be simulated by a spherically symmetric fluid flow. The 
Minkowski geometry in spherical coordinates with origin at the center of the fluid flow has the form 

r\ = -c 2 dt 2 + dr 2 + r 2 dfl 2 . (5) 

where as usual 

dn 2 = dd 2 + sm 2 9 d(j> 2 (6) 

is the 2-sphere metric. The fluid 4- vector can then be written in the form = (cosh — sinh (, 0, 0) where 
C = arctanh(t>/c) = ((r) is the boost parameter and v = \v\ = —v r is the velocity. We have chosen a minus sign 
for the radial component of it M to make the fluid flow inward directed. The Gordon metric becomes 

g = r] + fc(coshCcdt + sinh£dr) 2 . (7) 

It turns out that this metric will represent a Schwarzschild black hole with mass M if the fluid has a constant 
refractive index and its velocity is given by 

(8) 



/1 + (71 2 -1) 

'g 

where r g = 2GM/c 2 is the gravitational radius. The corresponding fluid 4- velocity has the components 



t > 2 -/ r /I"/ 



1 Since we are dealing with two different geometries, r\ and g, it is important to specify which metric is involved in raising and 
lowering of indices. 
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where / = 1 — r g /r. The velocity of light in the medium is c m = c/n. Comparing this with the fluid velocity 
itself, we find as expected that the fluid becomes superluminal at r = r s , i.e. at the black hole horizon. Of 
course the problem of the physical realization of the associated fluid flow remains to be addressed. The only 
possibility for an experimental realization of the horizon seems to be to have n large. In that case we can make 
the approximation v ~ c m y / r s /r. As it turns out, if the speed of sound of the fluid is constant, this flow also 
generates an acoustic black hole with a conformally Schwarzschild metric [Hj- To see what the flow given by 
equation JSJ is like we must solve the continuity equation for a stationary spherically symmetric fluid 

V • (pv) = ^-(r 2 pv) = . (10) 
dr 

This gives the radial dependence of the density as 

p=5 /l + (n 2 -l)-f , (11) 

where K is an integration constant. We see that the flow is compressible. 

We have defined a fluid flow which gives rise to an effective dielectric metric which coincides with the 
Schwarzschild vacuum geometry. To make this more clear we write down the solution in its explicit form (from 
now on setting c = G = 1) 



g = -rT 2 jdt 2 + 2n~V(l ~ f)( n2 ~ I) dtdr + I 1 + ™~ 2 (! - /)] dr 2 + r 2 dft 2 , (12) 

recalling that / = 1 — 2M / r and that n is constant. The usual form of the Schwarzschild metric can be recovered 
by the coordinate transformation 



dt = ndt s + rW(l-f)(n 2 -f)dr , (13) 
where tg is the usual Schwarzschild time coordinate. The metric then takes the familiar form 

g = -fdtl + /- 1 dr 2 +r 2 dQ 2 . (14) 
For reference we also give the integrated form of the transformation (|13(l 



t = nt s+ + 2Mnln ft - ^ ^ + n ^U) - f . (15) 

At infinity, where / = 1, the metric ()12J) reduces to 

g = -n- 2 dt 2 + dr 2 + r 2 dfl 2 . (16) 

This is the flat Minkowski space-time as it should be. The time component of the metric is rescaled accounting 
for the slower speed of light in the medium. 

We conclude this section by showing that these coordinates are related to both the Kerr-Schild and the 
Painleve-Gullstrand coordinates. First, we note that the metric H12(l is really a family of metrics parametrized 
by n. As it appears in l|12|l . the metric is defined for all values of n in the range n > 1, even though the 
underlying simulation is only defined for n > 1. Setting n = 1, the metric becomes 

9 = V + (dt + dr) 2 , 17 

r 

which is the Kerr-Schild form of the Schwarzschild geometry. This is very natural since the fluid velocity tends 
to the speed of light in the limit n 1 as can be seen from iJSJ. The metric also has a limit as n — > oo. To 
perform this limit we must first introduce a new time variable T — t/n in (|12|l . The metric then reads 



/ 2M\ , 9 2M ( 1 2A/\ , , / 2M\ 

= - 1 dT 2 + 2J 1 - + dTdr+ 1 + -5-1 dr 2 + r W . (18) 

\ r J y r \ n A n^r J \ n^r J 
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It is now possible to take the limit n->oo and the result is 




1 _ t£L ) d T 2 + 2\ — dTdr + dr 2 + rW . 




(19) 



This is just Schwarzschild's metric in Painleve-Gullstrand coordinates. Finally, introducing the parameter 
s = 1 — k = n~ 2 taking values in the interval < s < 1, the metric becomes 



This metric form interpolates between the Painleve-Gullstrand coordinates at s = and the Kerr-Schild coor- 
dinates at s = 1. 

3 Discussion 

In the metric given by l)12|l . the interpretation of the time coordinate t is that it is the time of the lab in which 
the simulation is taking place. This time is different from the Schwarzschild time tg. The radial variable r 
plays two roles. First, it is the Euclidean radial distance in the lab. But it is also the Schwarzschild radial 
variable in the simulated black hole. Another thing to note about the simulation is that there are actually two 
background geometries, one which is the background in the absence of the medium with metric (jSJl. This is the 
lab geometry. There is also the background in the medium at rest with metric (|16[1 . Both these backgrounds 
are fiat and they are related by the time transformation t = nT . We can interpret T as an "electromagnetic 
time" in the medium at rest. This means that a clock based on light travel time which shows the time t outside 
the medium would show the "slower" time T inside the medium. 

As discussed in the introduction, the Gordon form of the metric is interesting in its own right. Having found 
the flow corresponding to a Schwarzschild black hole, a natural question to ask, is if it would be possible to do 
the same thing with Kerr's rotating black hole. In a wider perspective one could ask which space-times admit 
a Gordon representation. The model presented in this letter has the property that the refractive index of the 
fluid is constant. However, models with a non-constant n are also possible. A physical requirement for such 
models to serve as simulations would be the existence of a constitutive relation n — n(p). In the spherically 
symmetric case this entails no restriction since there is only one non-trivial variable, and a constitutive relation 
is then always given parametrically by n — n{r) and p = p(r). In situations with more than one non-trivial 
variable, however, one would need to consider the Einstein equations coupled to the continuity equation. 

A more complete account of this work, including the method of finding the solution, will be given in a 
separate report. 
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g = -fdT 2 + 2y/(l - /)(1 - sf)dTdr + (1 + s - sf)dr 2 + r 2 dtt 2 . 
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